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ABSTRACT 

We present ISO-PHOT spectra of the regions 2.5-4.9 jum and 5.8-11.6 fim for 
a sample of 45 disk galaxies from the U.S. ISO Key Project on Normal Galaxies. 
The galaxies were selected to span the range in global properties of normal, star- 
forming disk galaxies in the local universe. The spectra can be decomposed into 
three spectral components: (1) continuum emission from stellar photospheres, 
which dominates the near-infrared (2.5-4.9 jttm; NIR) spectral region; (2) a weak 
NIR excess continuum, which has a color temperature of ~ 10 3 K, carries a 
luminosity of a few percent of the total far-infrared dust luminosity Lfir,; and 
most likely arises from the ISM; and (3) the well-known broad emission features 
at 6.2, 7.7, 8.6 and 11.3 /im, which are generally attributed to aromatic carbon 
particles. These aromatic features in emission (AFEs) dominate the mid-infrared 
(5.8-11.6 /mi; MIR) part of the spectrum, and resemble the so-called Type-A 
spectra observed in many non-stellar sources and the diffuse ISM in our own 
Galaxy. The few notable exceptions include NGC4418, where a dust continuum 
replaces the AFEs in MIR, and NGC 1569, where the AFEs are weak and the 
strongest emission feature is [S IV] 10.51 //m. 

The relative strengths of the AFEs vary by 15-25% among the galaxies. How- 
ever, little correlation is seen between these variations and either IRAS 60 //m-to- 
100 yum flux density ratio i?(60/100) or the far-infrared-to-blue luminosity ratio 
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Lfir/Lb, two widely used indicators of the current star-formation activity, sug- 
gesting that the observed variations are not a consequence of the radiation field 
differences among the galaxies. We demonstrate that the NIR excess continuum 
and AFE emission are correlated, suggesting that they are produced by similar 
mechanisms and similar (or the same) material. On the other hand, as the cur- 
rent star-formation activity increases, the overall strengths of the AFEs and the 
NIR excess continuum drop significantly with respect to that of the far-infrared 
emission from large dust grains. In particular, the summed luminosity of the 
AFEs falls from ~ 0.2L F tr for the most "IR-quiescent" galaxies to ~ 0.1 L F m 
for the most "IR-active" galaxies. This is likely a consequence of the preferential 
destruction in intense radiation fields of the small carriers responsible for the 
NIR/AFE emission. 

Subject headings: galaxies: ISM — infrared: galaxies — infrared: ISM: contin- 
uum — infrared: ISM: lines and bands — infrared: dust 



1. Introduction 

The spectroscopic properties of galaxies in the wavelength range 2.5-12 /zm are much 
less well known than in the optical regime, where stellar emission dominates, and in the far- 
infrared, where thermal continuum emission from cool interstellar dust dominates. Prior to 
the Infrared Space Observatory (ISO; Kessler et al. 1996) mission, spectroscopic information 
on this wavelength region was available for only a few of the brightest galaxies and galaxy 
nuclei, which were observed with ground-based telescopes through a few spectral windows 
(e.g., Roche et al. 1991), or with the IRAS Low- Resolution Spectrometer (e.g., Cohen & Volk 
1989). However, this relatively unexplored spectral region contains important signatures of 
interstellar dust particles, in particular the broad emission features at 3.3, 6.2, 7.7, 8.6, and 
11.3 fim (Gillett et al. 1973) which in the past were often called the "Unidentified Infrared" or 
UIR features. These features are now known to be ubiquitous in the ISM of our Galaxy: they 
are seen in H II and PDR regions (e.g., Cesarsky et al. 1996; Roelfsema et al. 1996; Verstraete 
et al. 1996, 2001); planetary nebulae and circumstellar regions (e.g., Beintema et al. 1996); 
reflection nebulae (e.g., Boulanger et al. 1996, Uchida, Sellgren, & Werner 1998); and diffuse 



1 Bascd on observations with ISO, an ESA project with instruments funded by ESA member states 
(especially the PI countries: France, Germany the Netherlands, and the United Kingdom) and with the 
participation of ISAS and NASA. 
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cirrus clouds (e.g., Mattila et al. 1996; Lemke et al. 1998; Chan et al. 2001). Taken together, 
they radiate a significant fraction of the total IR emission from these regions (e.g., Puget & 
Leger 1989). 

While the exact identity of the carriers of the UIR features is still unresolved, it is 
generally agreed that the features arise from vibrational modes of a carbon-based, aromatic 
material, so we will refer to them hereafter as aromatic features in emission, or "AFEs." 
Possible candidate materials range from Polycyclic Aromatic Hydrocarbon molecules (here- 
after PAHs; Leger & Puget 1984; Allamandola, Tielens, & Barker 1985, 1989; Puget & 
Leger 1989) to Hydrogenated Amorphous Carbon grains (Duley & Williams 1981, 1988). 
The current picture is that the particles responsible for the emission, whether present as 
free molecules or being attached to larger grains, are transiently heated by the absorption of 
single UV photons to T ~ 10 3 K. The carriers of the AFEs may also play an important role 
in regulating the physical conditions in the ISM by contributing significantly to photoelectric 
heating of the gas (Bakes & Tielens 1994; Helou et al. 2001; however, for an alternate point 
of view, see Chan et al. 2001). 

The spectral region also hosts a possible near-infrared, non-stellar continuum emission 
with a color temperature of ~ 10 3 K, which was first detected in Galactic reflection nebulae 
(Andriesse 1978; Sellgren, Werner, & Dinerstein 1983; Sellgren 1984). Recent observations 
with COBE and IRTS hint that this near-IR continuum may also be present on larger scales 
in our Galaxy (Bernard et al. 1994; Tanaka et al. 1996), but it has not been previously 
reported in external galaxies. 

With its unprecedented sensitivity and contiguous wavelength coverage, ISO made it 
possible to obtain 2.5-11.6 /zm spectra of large numbers of galaxies that were too faint to 
observe previously. As part of an ISO Key Project to study the physical properties of the 
ISM in galaxies (Helou et al. 1997; Dale et al. 2000), we obtained 2.5-11.6 /zm spectra using 
the PHT-S mode of ISO-PHOT (Lemke et al. 1996) for 45 galaxies. This survey covered 
the full range of morphological types of disk galaxies, SO to Im, that are powered by star 
formation. Similar ISO spectra of other types of galaxies can be found elsewhere: e.g., AGNs 
(Clavel et al. 2000), ultraluminous IR galaxies (Genzel et al. 1998; Rigopoulou et al. 1999), 
and galaxies of moderately low surface brightnesses such as the Magellanic Clouds (Reach 
et al. 2000; Sturm et al. 2000; Vermeij et al. 2002). The AFEs are seen in the spectra of 
most of these objects, but are absent or too weak to detect in elliptical galaxies (Lu & Hur 
2000; 2003; Athey et al. 2002) and extremely metal deficient dwarfs (e.g., Thuan, Sauvage, 
& Madden 1999). 

In Paper I (Helou et al. 2000), we presented spectra for a subsample of 7 galaxies, 
highlighted the detections of the AFEs and a near-infrared, non-stellar continuum emission, 
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and suggested that the averaged spectrum of a number of galaxies could provide a useful 
template for redshift determinations of distant star-forming galaxies. In the present paper, 
we explore more fully the diversity of the 2.5-12 /zm spectra of the galaxies in our sample. 
We present and compare the individual spectra of all the 45 observed galaxies, provide 
a quantitative analysis of various spectral components, evaluate variations from galaxy to 
galaxy and possible statistical correlations with galaxy properties at other wavelengths, 
and discuss possible reasons for these trends. Throughout this paper, we use NIR and 
MIR to refer to the 2.5-4.9 fxm and 5.8-11.6 /im spectral regions respectively, and AFE(3.3), 
AFE(6.2), AFE(7.7), AFE(8.6) and AFE(11.3) for the corresponding individual features. 



2. Galaxy Sample, Observations and Results 
2.1. The Sample 

The galaxies studied in this paper are a subset of a larger sample observed for the Key 
Project. This parent sample consists of 69 "normal" galaxies selected to capture the great 
diversity in the properties of galaxies in the local universe, especially in terms of the ratio 
of current to past star-formation rate. For each galaxy in this project we obtained at least 
one of the following: ISO-CAM images at 7 and 15 /im (see Dale et al. 2000), a sparsely 
sampled ISO-LWS spectrum between 43 and 200 /xm targeting fine-structure lines (Malhotra 
et al. 2001 and references therein), and a PHT-S 2.5-11.6 /xm spectrum as presented here. 

Table 1 lists the 45 Key Project galaxies for which we obtained a PHT-S spectrum. The 
position of each observation in the table was reconstructed from the ISO pointing history 
(IIPH) file, which in most cases was the same as the intended position. Nearly all of the 
observations were taken at the galaxy optical center, except for two off-center positions in 
NGC 1569 which correspond to emission peaks on the CAM images (Hunter et al. 2001). 
From Table 1 it can be seen that the blue luminosity L B ranges from 2 xlO 9 to 7 x 10 10 L &) 
the FIR-to-blue luminosity ratio Lfir/^b from 0.2 to 15, and the IRAS 60-to-100 /im color 
index i?(60/100) from 0.3 to 1.3. The heliocentric velocities [Col. (10)] were used to shift all 
spectra to a common (rest) frame for easier comparison. The median velocity of the sample 
is about 1800 km s -1 . 

Fig. 1 shows the distribution of the galaxies in a plot of logLpm/^B vs. i?(60/100). 
Most of the galaxies fall along a rough diagonal from the lower left to upper right corner of 
the diagram, a trend which we attribute to an increasing ratio of present-day star forming 
activity to the time-averaged star-formation rate in the past (e.g., Helou 1986). In this paper 
we assume that both axes in Fig. 1 are statistically valid measures of the global amount of 
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current star-formation activity in a galaxy, and thus of the UV-to-optical spectral shape 
of the radiation field. We define three galaxy subsamples: an "FIR-quiescent" subsample, 
represented by open squares in Fig. 1; an "FIR-active" group, represented by filled squares; 
and an intermediate subsample, shown as crosses. The few outliers in Fig. 1 include the 
compact galaxy NGC4418, which has higher i?(60/100) than 99% of the galaxies in the 
IRAS Bright Galaxy Sample (Soifer et al. 1989), and the low-metallicity irregular galaxy 
NGC 1569 which has low L FIR /L B . Most of the outliers also show peculiarities in their 
PHT-S spectra; we comment on these individually in §3.2. 



2.2. Observations and Data Reduction 

The PHT-S spectrometer has two 64-element linear Si:Ga detector arrays. The SS array 
covers the wavelength range 2.5-4.9 /im with spectral resolution element 0.04 /im while the 
SL array covers the interval 5.8-11.6 /im with a resolution element of 0.1 fim. The instrument 
has a 24" x 24" aperture on the sky, pointed with an accuracy < 2" (Kessler et al. 1996). For 
each galaxy, the SS and SL spectra were obtained simultaneously, through an aperture placed 
at the position given in Table 1. Sky reference observations were taken at symmetrically 
placed offsets of ± 150" from the galaxy center along a direction determined by the spacecraft 
roll angle at the time of the observation. Integration times were 512 seconds, split evenly 
between the galaxy and sky positions, with a duration of 64 sec per chopper step, except for 
the faint object IC860, for which we increased the total integration time to 2048 seconds. 

The spectra were derived from the Edited Raw Data using standard procedures in the 
PHOT Interactive Analysis package (PIA version 7; Gabriel et al. 1997), including deglitching 
at both ramp and signal levels, ramp slope fitting, signal averaging per chopper position, and 
sky subtraction of the average signal for the two sky reference positions. The flux calibration 
was performed using a signal-dependent "detector response function" obtained from chopped 
observations of calibration stars with known SED's. This "direct calibration," which was 
later incorporated into pipeline version 8.4, included an empirical correction for the signal 
loss due to a detector transient induced by the chopper switching between the source and 
the reference positions. 

Since we used point-source flux standards, our spectra represent the effective emis- 
sion corresponding to the integration of a normalized PHT-S beam profile over the surface 
brightness distribution of the source, divided by f ps f, the fraction of the point-source spread 
function within the PHT-S aperture (see Appendix A). For the SS detector pixels, the me- 
dian flux uncertainty a^ dian ~ 25 mJy and depends only weakly on the source flux. In 
contrast, u^ dian varies from 15 mJy to 45 mJy over the sample and scales roughly linearly 
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with source brightness. According to the data validation report from pipeline version 8.4, the 
absolute and relative flux uncertainties are on order of 10% for point sources. As described 
in Appendix A, we performed an independent check on the flux calibration by comparing 
our PHT-S data with the CAM imaging data of Dale et al. (2000); these agree to better 
than 18% for 5-8.5 /im and < 25% for the 4-5 /im region. 

2.3. Near-Infrared Photometry 

We used near-infrared images from the Large Galaxy Atlas of the 2-Micron All Sky 
Survey, 2MASS (Jarrett et al. 2003) to derive integrated J (1.25 /xm), H (1.65 jum), and Ks 
(2.17/im) magnitudes for each galaxy, appropriate to the placement of the PHT-S aperture 
on the sky during the ISO observations. Since there is not yet a finalized conversion formula 
from the 2MASS magnitude scale to flux densities, we normalized our 2MASS magnitudes 
to the multi-aperture observations of four early-type galaxies (NGC 4374, NGC 3379, NGC 
5866, and NGC 1326) published by Frogel et al. (1978) and Persson, Frogel, & Aaronson 
(1979), using the magnitude-to-flux conversions of Wilson, Schwartz, & Neugebauer (1972). 
To adjust our 2MASS photometry to the magnitude scale in these references, we found it 
necessary to add the following offsets: (0.028 ± 0.009) m at J, (0.052 ± 0.004) m for H, and 
(0.079 ± 0.010) m to convert from K s to K (2.2 ^m). The resulting 2 MASS fluxes are used 
to normalize our PHT-S spectra in §3 and §4 to unveil the NIR excess continuum in disk 
galaxies. This requires no significant zero-point offset between the 2MASS and PHT-S flux 
scales. We show in Appendix B that this is indeed the case. 

2.4. Results 

The sky-subtracted, rest-frame PHT-S spectra are presented in Figures 2a-2e. Table 2 
summarizes the PHT-S aperture coverage factor p (see Appendix A), the mean continuum 
flux density at 4 /im, and a mean flux density for each of the AFEs. These mean flux densities, 
defined in the footnotes to Table 2, are derived after shifting the spectrum to the rest frame 
and resampling the data at the PHT-S detector wavelengths by a linear interpolation between 
the two nearest data points. Note that the mean flux density for AFE(11.3) is basically taken 
over only the blue side of the feature, since part of the red side fell outside our wavelength 
coverage. The last 2 lines in Table 2 provide information on two elliptical galaxies, NGC 3379 
and NGC 4374, whose PHT-S spectra were discussed by Lu & Hur (2000; 2003), and which 
served as comparison objects for the disk galaxies (§4.1). 
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3. The Spectra 

3.1. General Characteristics 

The great majority of the observed galaxies show qualitatively similar PHT-S spectra. 
The MIR part of the spectra is typically dominated by the prominent AFE features at 
6.2, 7.7, 8.6, and 11.3 /im, which have relative strengths and profiles similar to those of 
"Type A" Galactic sources as described, for example, by Geballe (1997) and Tokunaga 
(1997). These features are superposed on a NIR continuum dominated by the emission from 
stellar photospheres. Of all the spectra, only that of NGC4418 display MIR structures that 
are qualitatively different from a Type A spectrum. We discuss in more detail this unique 
spectrum and a few others in §3.2. 

The similarity of most of the spectra in the range 5.8-11.6 /xm suggests that it is useful 
to derive an averaged spectrum to serve as a template for normal, star-forming galaxies. We 
derived such an average spectrum from 40 of the 45 galaxies in this study, omitting only the 
atypical galaxies NGC 4418 and NGC 1569 (see above, and §3.2), and the three galaxies with 
the lowest S/N ratios in their spectra (NGC 3705, NGC 4519, and NGC 7418). The average 
spectrum was obtained on a wavelength pixel-by-pixel basis from the resampled, rest-frame 
individual spectra (see §2.4) normalized at some fiducial wavelength. No correction for 
redshift was made to the JHK points described in §2.3, since the effects are negligible for 
these broad photometric bands. 

The averaged spectra resulting from two choices of the fiducial wavelength for normal- 
ization are shown in Fig. 3. The open squares depict the spectrum obtained by normalizing 
the individual spectra at J, while the thick solid curve corresponds to normalizing by the in- 
tegrated flux of AFE(7.7). In both cases we used 1/<t 2 - weighted averaging, where a is either 
^median or ^median d e p enc ii n g on which detector array segment the pixel under consideration 
belongs to; the JHK data points were averaged using the same a-weights as used for the 
PHT-SS data. 

The error bars in Fig. 3 represent s/y/n, where s is the standard deviation of the 
normalized fluxes and n is the number of galaxies included in the average (due to the redshifts 
of individual galaxies, n < 40 for some points near the array edges). For a given pixel, s is 
given by 

s 2 = (Zwiff- <f> 2 EwOA^Eti;*), (1) 

where fi and Wi are respectively the flux density and weight from the ith sample galaxy, and 
< / > is the average flux of the pixel. Therefore, these error bars reflect mainly the variation 
of spectral shape within the sample. Note that the two normalization methods have different 
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biases: normalization at 7.7 /im gives greater weight to individual spectra which are more 
dominated by emission from the ISM (e.g., the AFEs), whereas normalization at J gives 
greater weight to galaxies with a larger proportional contribution of starlight. 

The averaged spectra in Fig. 3 can be compared with the "template" spectrum presented 
in Fig. 2 of Paper I. The latter was a straight average [i.e., with Wi = 1 in eq. (1)] over a 
subset of 28 galaxies. In addition, a different PHOT calibration and a different normalization 
parameter (the strength of the 6.2 /iin AFE) was used. The main differences are limited to 
A < 3 /im; the MIR region of the template spectrum from Paper I is very similar to those 
derived here if they are all normalized in the same way. 

We also compare our average disk galaxy spectra with those of the two elliptical galaxies, 

NGC3379 and NGC4374. The PHT-S spectra of these two El galaxies, reduced from the 

ISO archive data by Lu & Hur (2003), are shown in Fig. 3 as a thin solid line and dotted 

line respectively. At a heliocentric redshift of Vh = 911km s _1 , NGC 3379 (M 105) is one 

i 

of the nearest normal giant elliptical galaxies, with a classical r4 profile (de Vaucouleurs & 
Capaccioli 1979). NGC 4374 (M84), at v h = 1060km s -1 , is known to host a central radio 
source (cf. Bridle & Perley 1984), and contains at least some interstellar matter, indicated 
by the presence of dust lanes and IRAS detections at 60 and 100 /im. Nevertheless, these 
two ellipticals have quite similar PHT-S spectra, nearly featureless and falling roughly as a 
Rayleigh- Jeans law from 1.25 to ~ 7/im. An apparent flattening beyond 8/xm may be due 
to circumstellar dust emission (Knapp, Gunn, & Wynn- Williams 1992; Athey et al. 2002). 

It can be seen from Fig. 3 that the spiral galaxies show clear excess emission compared 
to the ellipticals at all wavelengths > 2.2 /im. In Table 3, we list the numerical values of 
the two average disk-galaxy spectra in Fig. 3, as well as the average for the two E galaxies. 
These averaged spectra, because of their high S/N, offer a more sensitive way to study the 
profiles of the brighter AFEs and to look for weak features that might be buried in noise in 
individual spectra. 

Table 4 summarizes the features that are identifiable in the averaged spectra in Fig. 3 and 
were discussed in Paper I. The equivalent widths of some of weaker features were estimated on 
the AFE(7.7)-normalized spectrum and are given in the last column of Table 4. The 3.3 /im 
feature, which is not apparent in most of the individual galaxy spectra, is easily recognized in 
the averaged spectrum, and the narrow feature at around 4.03 /im is probably H I Bra. The 
hump near 7.0 /im could be either [Ar II] 6.99 /xm and/or the 6.9 /im dust feature discussed 
by Bregman et al. (1983) and Cohen et al. (1986), with possibly a small contribution from 
H 2 6.910/xm v—0-0 S(5). The small bump at around 10.6/xm could arise from [S IV], which 
is seen most distinctly in the post-starburst galaxy NGC 1569. Finally, the asymmetric 
appearance of AFE(7.7) is probably a result of the blending of two unresolved peaks at 7.6 
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and 7.8 /xm; such substructure has been seen in some Galactic sources (e.g., Bregman 1989; 
Roelfsema et al. 1996), and is partially resolved in some of the highest-S/N spectra in this 
study (e.g. NGC4194 and IC883, see Fig. 2). 

In order to assess whether the the MIR spectral shape varies with star-formation activity, 
we compared the average spectra for the FIR-quiescent and FIR-active subsamples defined 
in §2.1, which include 14 galaxies each. The average spectra for the J-band normalization 
described above, aligned at J for comparison, are shown in Fig. 4a as the solid curve (FIR- 
quiescent) and open squares (FIR-active). It can be seen from the figure that both the NIR 
continuum and AFEs for the FIR-active galaxies lie well above those of the FIR-quiescent 
subsample. On the other hand, when the averages are normalized by AFE(7.7) and aligned 
at 7.7 /xm (Fig. 46), the FIR-quiescent spectrum lies above the FIR-active one only in the 
NIR region. (The small difference between the two spectra in Fig. 46 in the region 9-10 /xm is 
significant only at the 1.5 a level.) Note that the apparent "dips" around 3.2, 3.7 and 4.1 /xm 
seen in Fig. 46 are probably not real, but are due to the AFE(7.7) normalization, which yields 
an overly strong weighting for the (noisy, and in some cases negative) short-wavelength data 
points of galaxies with strong AFEs and faint NIR continua. 

These plots demonstrate that the only significant spectral difference between the two 
subsamples lies in the NIR, and can be attributed to different contributions from starlight. 
In other words, the spectral shape of the MIR AFEs in galaxies remains largely independent 
of star-formation activity. 

3.2. Remarks on Individual Spectra 

3.2.1. NGC 520 

NGC520 (Arp 157), classified as an intermediate-stage merger by Hibbard & van 
Gorkom (1996), is as radio- and infrared-bright as NGC 4038/4029 (Arp 244; the Anten- 
nae). Numerical simulations suggest that NGC 520 is a merger remnant resulting from two 
disk galaxies which began colliding about 300 million years ago (Stanford & Balcells 1991). 
Our PHT-S aperture was placed on the main optical ridge, ~ 14" north of the emission peak 
in the 10/xm image of Bushouse, Telesco, & Werner (1998). This peak is presumably the 
dust-obscured nucleus of the brighter galaxy in the pair (see Fig. 6 of Stanford & Balcells 
1990). The PHT-S spectrum therefore included both this obscured nucleus and parts of the 
disks of both galaxies. 

The most striking aspect of the spectrum of NGC 520 is the weakness of AFE(8.6) and 
AFE(11.3) relative to AFE(7.7). On the other hand, AFE(6.2) has normal (Type A) strength 
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relative to AFE(7.7). In the context of the PAH model, both AFE(8.6) and AFE(11.3) arise 
primarily from C-H bending modes, while AFE(6.2) and AFE(7.7) come from C-C modes. 
Thus, a plausible conjecture is that the compression of the ISM due to the interaction of the 
merging galaxy disks might have led to a high degree of dehydrogenation of PAHs, resulting 
in relatively weak 8.6 and 11.3 /im features. 

3.2.2. NGC 1569 

NGC 1569 is a nearby Magellanic-type irregular galaxy with a metallicity ~ 30% of 
solar. It experienced a strong burst of star formation as recently as a few million years ago 
(Israel 1988; Israel & van Driel 1990; Greggio et al. 1998). We obtained PHT-S spectra at 
three positions on the galaxy disk. Position C (see Table 1) is on the nucleus, which also 
hosts one of the two superluminous young star clusters in this galaxy, Cluster B (Abies 1971; 
Arp & Sandage 1985). Positions NW and SE are located (1670 W, 5" N) and (16'.'6 E, 7" S) 
of the nucleus respectively. These correspond respectively to the NW and SE peaks on our 
CAM LW2 image (Hunter et al. 2001) as well as to H II regions No. 2 and No. 7 in Table 3 
of Waller (1991). 

All three spectra show a strong, unresolved emission line at 10.55 (±0.05) /mi, which 
presumably is [S IV] 10.51 /zm. We may also have detected [Ar III] 8.99 jum, in the NW 
and SE spectra. If present, it is a factor of 5-10 weaker than [S IV] 10.51 /im, indicating 
a relatively high effective temperature for the illuminating radiation field (e.g., Rank et 
al. 1978; Rubin 1985). This is consistent with the conclusions of Hunter et al. (2001), who 
infer T e ff = 40,000 K based on the Key Project LWS spectroscopy and other data. The 
large inferred number of early-type O stars indicates recent vigorous star-formation activity 
in NGC 1569 (Hunter et al. 2001). 

The unusual strength of the ionic lines relative to the dust emission in NGC 1569 
carries over to longer wavelengths. Its value of L([0 III] 88 //m)/LpiR, where L([0 III] is the 
luminosity in the [O III] 88 /im line emission, is higher by an order of magnitude or more 
(1.0-1.5 dex) than that of nearly all of the other galaxies in the Key Project sample (see 
Fig. 6a of Malhotra et al. 2001). The only other galaxy with such high-contrast ionic line 
emission is IC 4662, for which, unfortunately, we did not obtain a PHT-S spectrum. Recall 
also that NGC 1569 was an outlier in Fig. 1, having unusually low L FIR /L B for its IRAS 
color index 72(60/100). Whether this is due to the recent starburst or is a consequence of 
low metallicity is unclear, but the weak dust compared to the gas emission suggests a low 
dust-to-gas ratio in NGC 1569, and therefore possibly a high mean free path for far-UV 
photons as well as a hard UV radiation field. 
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3.2.3. NGC 44 18 

The PHT-S spectrum of NGC 4418 is markedly different from those of nearly all the 
other sample galaxies: the AFEs are not seen, and the MIR spectrum is dominated by an 
apparent emission plateau extending from 6-9 /im, which has been interpreted as arising 
from a broad continuum upon which is superposed a deep 10 /im silicate absorption feature 
(Roche et al. 1986). 

NGC 4418 was previously noticed to be an extreme object by Malhotra et al. (1997), 
who pointed out its extreme deficiency in the [C II] 157 /im line. It is a good example of 
the now well-established trend, that L([C II] 157/zm)/L F iR to Lfir decreases as the FIR 
luminosity and intensity of the dust-heating radiation field increases (Malhotra et al. 1997; 
Luhman et al. 1998). The nature of the illuminating radiation field in NGC 4418 is unclear, 
although its FIR color index, _R(60/100) = 1.3, is the highest ("warmest") in the full Key 
Project sample (see Table 5 of Malhotra et al. 2001). NGC 4418 is a very compact source 
at radio wavelengths, with a radius of < 0'/5 (~ 53/i _1 pc) at 20 cm (Condon et al. 1990), 
so our PHT-S spectrum enclosed most (84%) of its MIR emission. On the one hand, our 
PHT-S spectrum of NGC 4418 resembles the ISO spectrum of an ultra-compact H II region 
in M17 (Cesarsky et al. 1996), which would be consistent with the suggestion that NGC 4418 
contains an intense nuclear starburst of very high optical depth. On the other hand, Spoon et 
al. (2001) report seeing absorption features from ices in the PHT-S spectrum, and interpret 
NGC 4418 as a dust- and ice-enshrouded active galactic nucleus (AGN). 



3.2.4. IC 860 

IC 860 is also compact, with a 20 cm diameter of < 0V4 (or lOO/i^pc; Condon et al. 1990), 
consequently, our PHT-S aperture included about 81% of its MIR emission. This galaxy has 
been generally classified as a non-AGN (e.g., Leech et al. 1989). Given its compact size, warm 
FIR color of i?(60/100) = 0.94, and high value of logLpm/^B = 1.05, it is plausible that 
IC 860 harbors a nuclear starburst. While the AFEs in IC 860 have typical relative strengths, 
they have a low collective intensity relative to the FIR dust emission (see §4.3). IC 860 was 
also one of the extremely "[C II]-deficient" galaxies noted by Malhotra et al. (1997). If the 
carriers of the AFEs are major contributors of photoelectrons that heat the gas, these two 
observations taken together are consistent with the destruction of a substantial fraction of 
the AFE carriers by an intense UV radiation field. 
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3.2.5. NGC 5866 

NGC5866 is a nearby, early- type (SO), edge-on disk galaxy, with a dust lane visible in 
optical and the "coolest" FIR color index in the entire Key Project sample, i?(60/100) = 0.3 
The CAM LW2 image (Malhotra et al. 2000) shows that the emission at 7 fim arises mainly 
from the edge-on disk of the galaxy. Our PHT-S aperture was centered roughly on the 
nucleus, and enclosed about 60% of the total flux of the galaxy. 

The PHT-S spectrum of NGC 5866 shows a strong continuum rising towards shorter 
wavelengths, consistent with the picture that the global emission from this galaxy is strongly 
dominated by stellar photospheric emission. There appears to be a broad emission feature 
that peaks at 7.9 /iin, instead of the usual AFE(7.7/xm) feature. The 8.6 fxm feature is weak, 
if present, and the 11.3 /im feature is much wider than usual. Taken together, these AFEs 
are unusually weak relative to the FIR dust emission as compared to the other FIR-quiescent 
galaxies in the sample (see §4.3 below). 

4. Analysis and Discussion 

4.1. Continuum Emission at 3 to 5 Microns 

Spectral synthesis models indicate that as long as the luminosity of a galaxy is dom- 
inated by relatively old stellar populations (i.e., older than a few Gyrs), the shape of its 
stellar continuum spectrum in near-infrared remains largely independent of the details of its 
star formation history (e.g., Bruzual & Chariot 1993). This is supported by results from 
2MASS, which show that normal galaxies display a much smaller color dispersion in the 
near-infrared than in the optical (e.g., Jarrett 2000). We therefore assume that the averaged 
elliptical galaxy spectrum (§3.1) is representative of the intrinsic spectral shape of the stel- 
lar component in our disk galaxies. It is uncertain whether the circumstellar dust emission 
longward of 8 /xm can also be scaled in this way, but for disk galaxies, the emission is much 
fainter than the ISM emission in this spectral region. 

We first investigated whether internal reddening by dust can account for the observed 
NIR excess emission. In Fig. 5, the same 40 galaxies 7 included in the average spectrum (see 
§3.1) are shown in two color-color plots: logarithmic flux-density ratios K/J vs. H/J, and 



7 In §4 we use galaxies drawn from this 40-galaxy sample only. In comparing PHT-S fluxes with the FIR 
fluxes (i.e., in Tables 5 and 7, and Fig. 9), we further exclude those few without a PHT-S aperture coverage 
factor p in Table 2. 
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f (4 fxm) I J vs. H/ J, respectively. The galaxies are coded according to their degree of FIR- 
activity in the same way as in Fig. 1, and the elliptical galaxies are shown as open circles. We 
considered two geometries for the dust extinction: a foreground dust screen (shown as solid 
lines in Fig. 5); and the case of uniformly mixed stars and dust (shown as dotted curves). 
In both cases, we used a near-infrared reddening law of A(X)/A(J) = (\/1.25fim)~ 1J , after 
Mathis (1990). We assumed zero reddening for the ellipticals. The tick marks indicate Ay = 
1,2, and 3 m respectively, for the dust screen model; in the star/dust mixed case, one requires 
about 3 times as much total dust column density in order to produce the same amount of 
reddening as for the foreground screen. Apart from this scale factor, the two cases produce 
very similar trends in Fig. 5 within the parameter space occupied by our galaxies. As one can 
see from the figure, most of the galaxies lie above the reddening lines, particularly in the lower 
panel, f(4 /zm) / J vs. H/ J. Therefore, it appears that dust reddening cannot account for 
the color differences between our sample of disk galaxies and the reference elliptical galaxies. 
There must be an additional continuum component which becomes more prominent at the 
longer wavelengths. 

We set out to reconstruct the character of this "NIR excess continuum" as follows. We 
assumed that the J and H emission in the disk galaxies arises only from stars. Next, we 
subtracted the averaged elliptical galaxy spectrum in Table 3 from the data for each disk 
galaxy, after reddening it to the H/ J color of the latter using the foreground screen case. 
(For the few cases where the observed H/ J color is slightly bluer than that of the ellipticals, 
we take Ay = 0.) Finally, to increase the S/N ratio, we averaged the residual spectra for the 
entire galaxy sample. The result is shown in Fig. 6a. Also shown is a modified black-body 
curve for T = 750 K and a A~ 2 emissivity law, which fits the general shape of the curve quite 
well. The emission is clearly much broader and stronger than AFE(3.3 /im), and therefore is 
unlikely to be simply the wings of this feature. 

The residual spectrum shown in Fig. 6a represents a lower limit to the typical intensity 
of this NIR excess continuum, because of our assumption that it does not contribute at H. 
In order to set an upper limit to such a contribution, we considered the effects of assuming 
that only the J-band flux has no contribution from non-stellar sources, and derived the 
residual spectrum that results from subtracting the elliptical galaxy spectrum uncorrected 
for different reddening between the disk and E galaxies. This yields the average residual 
spectrum shown in Fig. 66, which is well fit by a modified black-body curve with T = 10 3 K 
and a A -2 emissivity, only slightly hotter than the curve in Fig. 6a. Both approaches therefore 
indicate that the color temperature of the non-stellar NIR continuum is close to 10 3 K. The 
fractional energy contained in the NIR spectral region is about 10% in Fig. 6a and about 
17% in Fig. 66. 
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Next we consider whether the NIR excess continuum is correlated with AFE emission. 
In Fig. 7, we plot the ratio of the average 4 fim emission (taken to be representative of the 
NIR excess emission) to the mean flux density of AFE(7.7), against logLpiR/^B- (The two 
methods of removing the stellar continuum, discussed above, lead to more or less the same 
4/zm fluxes.) The strength of AFE(7.7/im) relative to the J-band flux increases by nearly 
an order of magnitude within the sample (see Table 2), as Lfir/£b increases. On the other 
hand, there is no net trend in Fig. 7, indicating that the NIR excess correlates strongly with 
the MIR emission for galaxies with a wide range in dust content and luminosity. We view 
this as strong evidence that the NIR excess continuum originates in the ISM, as opposed to 
being (scattered) starlight from late-type stars or circumstellar dust emission. However, the 
scatter in Fig. 7 is greater than can be accounted for by measurement errors, suggesting that 
there are some real variations among galaxies in the ratio of the NIR excess to AFE(7.7). 
The most deviant point is the early-type galaxy NGC 5866, which, as mentioned in §3.2.5, 
has very weak AFE features. 

This non-stellar dust continuum is probably the same component that has been detected 
in reflection nebulae and the large-scale ISM of our own Galaxy (see §1). In fact, Sellgren 
et al. (1985) showed that, for the reflection nebulae NGC 7023 and NGC 2023, the 4 /im 
emission surface brightness is between 1/6 to 1/10 of that at the peak of AFE(7.7). This 
range is quantitatively consistent with our results in Fig. 6. It has been suggested that such a 
NIR continuum might be due to electronic fluorescence or a quasi-continuum of overlapping 
bands of PAH molecules (e.g., Allamandola, Tielens, & Barker 1989). While further work 
is needed in order to determine whether the AFEs and NIR continuum emission arise from 
exactly the same macromolecules or material, it seems clear that the carriers of these two 
components are closely related. 

In Table 5 we summarize the luminosity in the NIR region relative to the FIR dust 
emission, for each of the three subsamples defined in §2.1. The total NIR luminosity is 
5-18% Lfir and arises mostly from stellar photospheres. If we correct for the stellar contri- 
bution using the two methods described above, the non-stellar NIR continuum contributes 
~ 3 — 4% Lfir for the most quiescent galaxies in our sample, and 1 — 2% Lfir for the most 
FIR-active ones. These estimates are on the same order as the COBE results for our own 
Galaxy (Bernard et al. 1994). 

4.2. Variations in the Mid-Infrared AFEs 

Some variations have been observed among Galactic sources in the profiles of the AFEs 
(e.g., Roelfsema et al. 1996; Peeters et al. 2002) and the ratios of one feature relative to 
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another (e.g., Joblin et al. 1996; Lu 1998; Vermeij et al. 2002), but the physical implications 
are still not fully understood. In particular, it appears that relative strengths of the fea- 
tures are quite insensitive to the intensity and color temperature of the local radiation field 
(e.g., Uchida et al. 2000; Chan et al. 2001). Likewise, while we find similar variations in the 
AFEs in our sample of galaxies, we have not found any statistically significant correlation 
between these variations and Lfir/^b, #(60/100), optical morphology, the 7 //m-to-15 /im 
flux density ratio (Dale et al. 2000), or the mean surface brightness at 7 or 15 /im within an 
isophotal ellipse containing 50% of the infrared flux. 

In Fig. 8 we plot the ratios of the strengths of the 6.2, 8.6, and 11.3 /zm AFEs to 
AFE(7.7) against Lfir/I/b for the 40 galaxies, which all show clearly discerned AFEs. We 
hereafter abbreviate the mean F\ of a feature by its wavelength, e.g. (6.2) for the mean flux of 
AFE(6.2). There is some indication of a slight decreasing trend in these plots. However, these 
trends are weak and probably not statistically significant. We computed the median feature 
strength ratios separately for the FIR-quiescent, intermediate, and FIR-active subsamples 
defined in §2.1. The values are given in Table 6 and plotted as large crosses in Fig. 8. It 
can be seen that the median values decline by less than 1 a s from FIR-quiescent to FIR- 
active, where a s is the r.m.s. dispersion of the entire sample (see the last row of Table 6). 
This differs from the conclusion of Lu et al. (1999) that the ratio (11.3) / (7.7) decreases with 
increasing L F i R / L B - That earlier study was based on an older PHT-S calibration, but more 
importantly, the stellar continuum was removed in a more simplistic way, by subtracting a 
power-law fit to the spectral intervals 3.5-5 fim and 9.5-10.5 /xm. That procedure had the 
effect of undersubtracting the continuum underlying the 11.3 //m feature in galaxies with a 
strong stellar component, thus producing an apparent trend. 

The scatter in Fig. 8 is significantly greater than the statistical error bars, implying that 
there are intrinsic galaxy-to-galaxy variations in the relative feature strengths. According to 
Table 6, the spread in the AFE ratios is ~ 15% for (6.2)/(7.7) and (8.6)/(7.7), and ~ 25% 
for (11.3) / (7.7). Within the framework of the PAH hypothesis, variations in the relative 
feature strengths have been attributed to the presence of different fractions of PAHs that are 
ionized and/or hydrogenated, as a consequence of differences in the ambient UV radiation 
field (e.g., Jourdain de Muizon et al. 1990; de Frees et al. 1993; Schutte et al. 1993; Langhoff 
1996). However, our data suggest that the situation may be more complex than a simple 
dependence on the UV radiation field. It is plausible that chemical processing may also play 
an important role in determining the local abundance and emission properties of the PAHs 
(e.g., Boulanger et al. 1990). 

The spectra in Fig. 2 also show some galaxy-to-galaxy variations in the AFE profiles. 
In Appendix C, we quantify one aspect of these variations by defining a "logarithmic slope" 
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S that measures the steepness of the feature profile on the short- or long-wavelength side. 
As we show in Table CI, the galaxy-to-galaxy variations in S range from 15-30% for all the 
AFEs, except for the short-wavelength side of AFE(8.6) which shows possibly much greater 
variations relative to the median value. This could either be due to intrinsic variations in the 
profile of AFE(8.6), or to a contribution from a weak feature around 8.2 /xm that varies from 
galaxy to galaxy. Such a shorter-wavelength feature was invoked by Verstraete et al. (2001) 
to fit high-S/N ISO-SWS spectra of Galactic sources, and an emission feature at ~ 8.2 /xm 
has also been observed in some post-AGB stars (Peeters et al. 2002). 

4.3. The Destruction of the AFE Carriers 

It was previously known that the AFEs are depressed in H II regions (e.g., Cesarsky 
et al. 1996). This has been widely attributed to the preferential destruction of the small 
AFE carriers relative to larger dust grains in these regions. IRAS data also indicate an AFE 
depression over galaxy scales, with the AFEs being more severely depressed in galaxies with 
warmer #(60/100) colors (Helou, Ryter, & Soifer 1991). 

Fig. 9 shows how Lafes, the summed luminosity of the four MIR AFEs (6.2, 7.7, 8.6, 
and 11.3 /mi), correlates with Lfir- The few galaxies without a PHT-S aperture coverage 
factor p in Table 2 are not plotted here. We indicate galaxies with values of p < 35% (for 
which the aperture correction is large and therefore less certain) by filled symbols, in order 
to distinguish them from the others (open symbols). For some of these galaxies, Malhotra 
et al. (2001, Table CI) have derived Go, the intensity of the far-UV interstellar radiation 
field in units of the radiation field in the solar neighborhood, based on a PDR model that 
depends on the observed fluxes in [C II] 158 /mi and [O I] 63 /xmand the flux of the total 
infrared emission. These galaxies are plotted as squares of three different sizes, where the 
smallest corresponds to 2.3 < logG < 3.4, the intermediate size to 3.5 < logG < 3.9, and 
the largest to 4.0 < logG < 4.8. 

For comparison, we also show in Fig. 9 the IRAS results of Werner, Gautier, & Cawlfield 
(1994) for multiple spatial positions in two Galactic H II regions, the Rosette and California 
nebulae (indicated by a dashed line), and the reflection nebula in the Pleiades (indicated 
by a solid line). These lines actually represent the logarithmic ratio of the IRAS 12 /xm flux 
(assumed to be proportional to our integrated AFE flux) to the FIR flux, and are shifted 
vertically in order to bracket most of the galaxy data points. The average trend of our 
galaxy points is steeper than that characteristic of the reflection nebula, but shallower than 
the trend for the H II regions. This comparison suggests that both H II regions and reflection 
nebulae contribute to the total AFE and FIR emission of galaxies. 
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In Table 7 we list median ratios of several different measures of the AFE overall intensity 
to the FIR flux, for the three subsamples. It can be seen from the table that the depression 
of AFEs relative to the FIR emission becomes significant only for the most FIR-active 
subsample; for the latter galaxies, Lafes/^fir has a value only 60% that of the other two 
subsamples. Furthermore, from Fig. 9 it can be seen that this AFE depression is more related 
to the hardness than the intensity of the radiation field: at a given 72(60/100), galaxies with 
higher values of G tend to lie closer to the dashed line. Since the flux ratio of [C II] to 
the AFEs has very little dispersion in this sample of galaxies (Helou et al. 2001), one does 
expect such a linkage between the diagnostics based on the fine-structure lines and the AFEs. 
However, this is unlikely just a mathematical linkage, for the two reference lines in Fig. 9 
clearly suggest that the (UV-hard) HII regions have a lower Z/afes/^fir than the (UV-soft) 
reflection nebulae at any given iZ(60/100). 

If the photodestruction process that weakens the AFEs in Galactic H II regions is 
also responsible for depressing the AFE emission on galaxy scales, our results imply that the 
mass spectrum of the dust particles could differ significantly between quiescent and starburst 
galaxies. 

As shown in §3.1, the MIR emission from disk galaxies is dominated by strong AFE 
emission over a wide range in such physical characteristics as the intensity of the interstellar 
UV radiation field, FIR luminosity, and temperature of the FIR-emitting dust. The few 
exceptions we find may be indicative of the extreme conditions required for the AFEs to fade 
from view. Representing the high-radiation-field limit, we have objects such as NGC 4418 
(which may be an AGN), IC 860, and NGC 1569 (which, if plotted in Fig. 9 despite its 
small value for p, would lie near the position of IC 860). For these objects, it is reasonable 
to explain the weakness of the AFEs as being the consequence of destruction of the AFE 
carriers by hard UV photons (or shocks), by analogy with the Galactic sources discussed 
above. 

It is less clear what causes the deficiency in AFE emission in the FIR-quiescent galaxy 
NGC 5866. This is the earliest type galaxy in Key Project sample for which we have a 
PHT-S spectrum; unfortunately, we did not obtain PHT-S observations for the other three 
E/S0 galaxies in our program (Malhotra et al. 2000). The lack of strong AFE emission 
in NGC 5866 could be due either to an actual deficiency of the AFE carriers or to a lack 
of energetic UV photons capable of exciting the features, or both. Note that our LWS 
observations of the [C II] 158 /zm and [O I] 63 fim lines indicate a particularly weak and soft 
UV radiation field in this galaxy (Malhotra et al. 2000). 
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5. Summary 

We present new ISO PHT-S spectra (2.5-11.6 jum), extended to 1.25 fim using JHK data 
from the 2MASS survey, for a sample of 45 disk galaxies that span the typical range in global 
properties of galaxies energetically dominated by star formation. PHT-S aperture coverage 
factors are also provided for most of the galaxies. We decompose the spectra into three 
constituents: (1) stellar continuum emission, which dominates at the shortest wavelengths 
(2.5-4.9 /iin; NIR); (2) a weaker and redder NIR "excess continuum"; and (3) the well-known 
aromatic dust emission features (AFEs) at 6.2, 7.7, 8.6, and 11.3 /im which dominate the 
long-wavelength half of the PHT-S spectra (5.8-11.6; MIR). Most of the galaxy spectra 
appear similar to each other, especially in the MIR region, and to an averaged "template" 
spectrum shown in Fig. 3. The most striking exceptions are NGC4418, which shows no 
AFEs but instead has a broad mid-infrared continuum upon which is superposed a strong 
10 /im silicate absorption feature; and NGC 1569, which has a very weak infrared continuum 
but an unusually strong [S IV] 10.5 /im line from ionized gas. 

The non-stellar NIR excess continuum has an average color temperature of ~ 10 3 K and 
a luminosity of a few percent of Lfir, where Lfir is the far-infrared (40-120 /im) luminosity 
attributed to large dust grains. This NIR continuum scales more or less linearly with the 
strength of the AFEs, suggesting that the NIR excess continuum originates in the ISM of 
galaxies (and not, for example, in late-type stars or circumstellar dust); and that the AFEs 
and NIR excess continuum arise from similar carriers. 

The profiles and relative strengths of the AFEs in the disk galaxies match those of "Type 
A" mid-infrared spectra, by far the most predominant type of mid-infrared dust emission 
pattern in the ISM of our own Galaxy. This resemblance suggests that the carriers of the 
Type-A AFEs are also prevalent in the ISM of other galaxies. The combined luminosity 
of the AFEs in the region 5.8-11.3 /im is 10-20% Lfir. The relative strengths of the AFEs 
vary on average by 15-25%. These observed variations, however, do not correlate with IRAS 
60 /im-to-100 /im flux density ratio, _R(60/100), or the far-infrared-to-blue luminosity ratio 
L FIR /L B , two commonly used indicators of global star-formation activity in galaxies. We 
interpret this as indicating that other factors, in addition to the present-day radiation field, 
affect the strengths and shapes of the AFEs. 

The ratios of both the AFEs and non-stellar NIR continuum to the FIR flux decrease 
systematically from the most quiescent galaxies to the most actively star-forming galaxies 
in the sample. We show that this is more related to the hardness than the intensity of the 
heating radiation field, likely a result of the AFE/NIR carriers being preferentially destroyed, 
relative to the larger dust grains responsible for the FIR emission, in regions of active star 
formation and intense UV radiation fields. This implies that the mass and size distribution 
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functions of the interstellar dust particles vary in environments with different levels of star 
formation activity. 
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ISO-CAM LW2 images are available for 42 galaxies in Table 1 from Dale et al. (2000), 
and for one additional galaxy, NGC 5866, in the ISO archive from the observations described 
by Vigroux et al. (1999). For each of these galaxies, we integrated the LW2 filter curve over 
the PHT-S spectrum to derive f p (LW2), a CAM LW2-equivalent flux density at 6.7/zm. 
The LW2 bandpass extends to about 5 fim on the short-wavelength side. For its wavelength 
coverage beyond the blue end of the SL array, we used a simple linear interpolation between 
the red end of the SS array and the blue end of the SL. 

We then centered the PHT-S aperture on the CAM image at the position of our PHT-S 
observation. Because a CAM filter- wheel jitter can affect the effective pointing of CAM 
images, a small residual positional uncertainty remains. After subtracting a constant sky 
level from the CAM image, we derived the quantity f c (LW2) by integrating the source 
surface brightness over the PHT-S aperture. As a result, we have 
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A. Flux Comparison of the PHT-S and ISO-CAM Data 
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SB(x, y)dxdy, 
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where SB(x,y) is the source surface brightness distribution function. We can express 
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f p (LW2) in a similar way: 

f p (LW2) = (l/f psf ) J B(x, y) SB(x, y)dxdy, (A2) 

where f ps f is the fractional point-source spread function for PHT-S and B(x,y) is the PHT- 
S beam profile (which has a value of unity near the aperture center). We have implicitly 
assumed that f ps f, B(x,y), and SB(x,y) are all effective values over the PHT-S detector 
pixels included within the CAM LW2 bandpass. The integration in eq. (Al) is over the PHT- 
S aperture only, while that in eq. (A2) is over the entire source. For a source that is sufficiently 
compact that it can be treated as a point source with respect to the PHT-S beam profile, one 
has f c (LW2)/f p (LW2) w f psf . For the PHT-S detector pixels relevant to the CAM LW2 
filter, f ps f a: 0.92. Therefore, for point sources we should have f c (LW2) / f p (LW2) ps 0.92, if 
there were no systematic flux offset between the two instruments. 

Fig. 10 is a plot of f c (LW2)/f p (LW2) as a function of p, the fractional CAM LW2 flux 
of a galaxy that falls within the PHT-S aperture, which is given by the ratio of f c (LW2) 
to the total LW2 flux of the galaxy. Notice that the scatter is larger for more extended 
sources (i.e., smaller values of p); this effect may be partly due to the residual positional 
uncertainty described above. In addition, f c (LW2) / f p (LW2) increases on average as p drops, 
presumably due to the beam size effect in eq. (Al). The average flux ratio for the 20 galaxies 
with p > 0.6 is (0.76 ±0.03) and is indicated by the dotted line in Fig. 10. This implies 
a possible systematic flux difference of no greater than 18% between the CAM LW2 and 
PHT-S results. 

A similar comparison was performed with the CAM LW1 (4-5 /im) data (Dale et 
al. 2000) on the 10 galaxies for which these data were available. The result is shown in 
Fig. 11. The scatter in this diagram is significantly greater than in Fig. 10, largely because 
the galaxies are much fainter in LW1. Nevertheless, a quantitatively similar result is found 
as for the LW2 comparison: the median flux ratio in Fig. 11 is 0.7, which implies a difference 
of about 25% in flux scale after taking f ps f into account. 



B. Zero-Point Difference between the PHT-S and 2MASS Fluxes 

Using the 2MASS K-band flux density and the average PHT-S 4/iin flux density as 
defined in Table 2, we show here that the zero-point offset between the 2MASS and PHT-S 
flux scales is insignificant. In Fig. 12a we plot f u (4 /jm) against f v {K) for our sample galaxies. 
The dotted line is a least-squares fit to all the data points, by minimizing their distances 
perpendicular to the fit, and intercepts the vertical axis at about (0.006 ± 0.005) Jy. However, 
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we believe that this overestimates the true zero-point offset. In fact, f u (A/j,m)/f v (K) should 
rise somewhat as logL FIR /L B increases, because the contribution of the non-stellar 3-5 /im 
continuum emission (see §4.1) has a greater influence on f v (A jum) than on f u (K). Indeed, 
Fig. 126 shows that there are more FIR-active galaxies at smaller K fluxes. A better estimate 
of the zero-point offset in flux is therefore given by the solid line in Fig. 12a, which is a similar 
least-squares fit but includes only those galaxies with log Lpi^/ Lb < 0.1 (which are shown 
as solid squares in both Figs. 12a and 126). This line has a vertical interception of only 
(0.002 ± 0.007) Jy. We therefore conclude that the zero-point offset between the PHT-S and 
SMASS flux scale is insignificant. 

C. Slopes of Feature Profiles 

For each AFE, we define a "logarithmic slope" parameter S that measures the steepness 
of the feature profile on the short or long- wavelength side, 

s='-f-±, (ci) 

I + 1 1 

where Iq is the flux density of the detector pixel nearest the peak of the emission feature in 
the rest-frame, and l\ is the flux density of a detector pixel down one side of the feature. 
For Iq we used the pixel corresponding to 6.216 /im in Table 3 for AFE(6.2), 7.616 /zm for 
AFE(7.7), 8.540 /im for AFE(8.6) and 11.263 for AFE(11.3). For simplicity we chose h 
at one of the detector pixels that was used to define the wavelength range of the feature (see 
Table 2). Two slopes can be defined for each AFE, £_ for the short- wavelength side and S + 
for the long-wavelength side, except for AFE(11.3), for which only S_(ll.3) can be defined. 

In Table CI, we list for each AFE: the median value for the 40-galaxy sample defined 
in §3.1 [Col. (2)]; the r.m.s. dispersion around this median [Col. (3)]; and a m , the median 
of the statistical measurement errors [Col. (4)]. While the ratio of Cols. (3) to (2) tells us 
about the variation in the slope, the ratio of Cols. (3) to (4) can be used to gauge whether 
this is significant. The ratio of Cols. (3) to (2) ranges from ~ 15-30% for both slopes of 
AFE(6.2) and AFE(7.7) and for S+(8.6), to 40% for S_(11.3), to 130% for 5_(8.6). The 
ratios of Cols. (3) to (4) show that, of these, only the variations in S-(8.6) are significant at 
3<7 m . Possible reasons for the observed variations in S-(8.6) are discussed §4.2. 
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Table 1. Galaxy Properties 
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Table 1 — Continued 



Galaxy 
(1) 


R.A. a 
(2) 


Dcc a 
(3) 


ROLL 6 

(4) 


ISO-TDT c 

(5) 


Morphology d 
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log L% 
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log L Fm /L B f 
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SA(s)a: pec 


10.09 


1.12 


0.76 


5541 



a J2000 Right Ascension and Declination in "hhmmss.s" and "ddmmss," respectively. 

b The ISO roll angle in degrees, measured in a counterclockwise manner from the celestial north to the spacecraft Z-axis which is 
perpendicular to the direction in which the chopper was operated for sky reference positions. This roll angle gives the sky orientation 
of the square aperture of PHT-S which is aligned with the spacecraft axes. 

c The TDT number of the ISO observation as it appears in the ISO archive. 

d Optical morphology taken from Dale et al. (2000) if available, otherwise from the RC3 catalog (de Vaucouleurs et al. 1991). 
e Logarithmic optical blue luminosity in solar units, derived as in Dale et al. (2000). 
^Logarithmic FIR-to-bluc luminosity ratio, derived as in Dale et al. (2000). 

g IRAS 60 ^tm-to-100 fim flux density ratio, derived from IRAS addscan fluxes [see Table 4 of Dale et al. (2000)]. 
h Hcliocentric velocity in units of km s — 1 . 
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Table 2. Mean Flux Densities F x around Selected Wavelengths" 
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u.oy i^u.uo ) 


NGC 3620 


0.75 


30.99 


27.11 


14.64 


2.33(0.11) 


3.95(0.04) 


6.38(0.03) 


2.89(0.02) 


1.78(0.05) 


NGC 3683 


0.47 


18.21 


14.35 


7.07 


1.25(0.11) 


2.09(0.04) 


2.90(0.02) 


1.33(0.02) 


0.91(0.06) 


NGC 3705 


0.20 


20.28 


14.60 


6.57 


0.85(0.09) 


0.50(0.05) 


0.39(0.03) 


0.20(0.02) 


0.07(0.05) 


NGC 3885 


0.82 


26.73 


19.72 


9.51 


1.83(0.14) 


1.87(0.03) 


2.57(0.03) 


1.30(0.03) 


1.40(0.04) 


NGC 3949 


0.24 


14.47 


10.04 


4.35 


0.47(0.10) 


0.77(0.04) 


1.06(0.02) 


0.51(0.02) 


0.55(0.05) 


NGC 4027 




7.83 


5.32 


2.30 


0.37(0.08) 


0.69(0.03) 


0.89(0.03) 


0.39(0.02) 


0.31(0.04) 


NGC 4102 


0.62 


45.15 


34.18 


17.04 


3.39(0.14) 


4.44(0.03) 


6.78(0.03) 


3.36(0.03) 


2.94(0.06) 


NGC 4194 


0.81 


11.84 


8.62 


4.27 


0.42(0.13) 


2.47(0.03) 


3.68(0.03) 


1.68(0.02) 


1.02(0.08) 


NGC 4418 


0.84 


6.68 


4.57 


2.09 


0.47(0.11) 


1.35(0.04) 


2.86(0.03) 


1.23(0.02) 


0.29(0.06) 


Mpp A A on 


U.Uo 


1 3 fi7 


o.DD 


3 AO 
O.DZ 


U.Oo^U. LZ ) 


fl AC\lf\ H3\ 
U.4U^U.UO ) 


n A 1 (f\ H3\ 
U.^tl^U.Uo ) 


n ooi'n n9\ 
\J.ZZ\\J.uZ ) 


n 33/0 n.Q"\ 
U.oo^U.Uo ) 


NGC 4519 


0.27 


4.62 


3.17 


1.31 


0.39(0.11) 


0.36(0.03) 


0.46(0.03) 


0.25(0.03) 


0.11(0.04) 


NGC 4691 




13.09 


8.70 


3.94 


0.91(0.16) 


1.94(0.07) 


2.69(0.04) 


1.39(0.04) 


1.38(0.06) 


IC 3908 


0.52 


11.51 


8.80 


4.42 


0.84(0.11) 


1.51(0.03) 


2.14(0.01) 


0.77(0.12) 


0.95(0.05) 


IC 860 


0.89 


5.30 


3.76 


1.73 


0.28(0.03) 


0.20(0.01) 


0.31(0.01) 


0.12(0.01) 


0.16(0.02) 


IC 883 


0.90 


3.88 


2.83 


1.53 


0.56(0.12) 


1.26(0.03) 


1.90(0.03) 


0.76(0.02) 


0.43(0.04) 


NGC 5433 


0.71 


7.62 


5.90 


2.99 


0.72(0.12) 


1.23(0.05) 


1.90(0.02) 


0.87(0.04) 


0.74(0.07) 


NGC 5713 


0.45 


16.58 


11.34 


5.26 


1.66(0.13) 


2.30(0.04) 


3.34(0.03) 


1.64(0.04) 


1.73(0.07) 


NGC 5786 


0.30 


11.42 


8.37 


3.80 


0.83(0.11) 


0.66(0.03) 


0.84(0.02) 


0.46(0.03) 


0.39(0.05) 


NGC 5866 


0.61 


64.78 


46.92 


22.32 


3.32(0.11) 


0.87(0.04) 


0.68(0.02) 


0.38(0.02) 


0.41(0.04) 


NGC 5962 


0.30 


22.20 


14.70 


6.71 


0.53(0.13) 


0.72(0.04) 


1.05(0.02) 


0.61(0.02) 


0.40(0.07) 


IC 4595 


0.34 


9.90 


8.58 


4.32 


0.75(0.11) 


1.04(0.04) 


1.45(0.02) 


0.68(0.02) 


0.70(0.05) 


NGC 6286 


0.69 


7.71 


6.41 


3.55 


0.90(0.14) 


1.56(0.06) 


2.44(0.02) 


1.00(0.02) 


0.87(0.06) 


NGC 6753 


0.33 


38.53 


28.44 


13.21 


1.64(0.10) 


1.19(0.02) 


1.40(0.02) 


0.70(0.02) 


0.82(0.03) 


NGC 7218 


0.36 


8.23 


5.58 


2.45 


0.41(0.10) 


0.58(0.03) 


0.78(0.02) 


0.35(0.02) 


0.35(0.04) 


NGC 7418 


0.12 


5.45 


3.70 


1.69 


-0.11(0.14) 


0.30(0.03) 


0.32(0.02) 


0.13(0.02) 


0.09(0.04) 
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Table 2 — Continued 



Galaxy 


p" 


J c 


H c 


K c 


4fj,m d 


(6.2) e 


(7.7) e 


(8.6) e 


(11.3) e 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


IC 5325 


0.22 


9.98 


6.84 


3.02 


0.23(0.08) 


0.48(0.02) 


0.59(0.02) 


0.24(0.02) 


0.23(0.05) 


NGC 7771 


0.52 


17.35 


13.41 


6.70 


1.16(0.10) 


1.80(0.04) 


2.74(0.03) 


1.34(0.03) 


1.13(0.07) 


Mrk 331 


0.76 


8.98 


6.72 


3.46 


0.80(0.12) 


1.82(0.04) 


2.73(0.02) 


1.32(0.03) 


1.22(0.07) 


NGC 3379 




119.31 


79.66 


34.94 


4.81(0.06) 


0.86(0.02) 


0.50(0.01) 


0.48(0.01) 


0.39(0.03) 


NGC 4374 




99.70 


71.52 


31.67 


3.77(0.06) 


0.61(0.01) 


0.41(0.01) 


0.40(0.01) 


0.36(0.02) 



a All flux densities arc in units of 10 -14 Wm — 2 (im -1 . To convert between this unit and Jy, we use 1.25, 1.65, 
2.2, 4/im as the effective wavelengths for Columns (3), (4), (5) and (6), respectively. The parenthesized value 
following a flux density is the statistical error of that flux. 

b PHT-S aperture coverage factor p as defined in Appendix A. 

C JHK flux densities within the PHT-S aperture. 

d The mean 4 (im flux density derived from a simple average of the flux densities in Jy of the 27 pixel channels 
between 3.4 and 4.4 fim. 

e The mean flux density of an AFE, defined as the integrated flux between Ai and A2 divided by (A2 — Ai), 
where Ai and A 2 arc respectively 5.98 ^m and 6.64 /nn for AFE(6.2), 7.20 /im and 8.22 ^mfor AFE(7.7), 8.22 (im 
and 9.23 ^mfor AFE(8.6), and 10.86 /um and 11.40 /zmfor AFE(11.3). 



Table 3. Rest-Frame Average Spectra" 



\(fim) 


Norm, at J 6 


Norm, by (7.7) c 


n d 


Ellipticals 6 


(1) 


(2) 


(3) 


(4) 


(5) 


1.250 


1.000(0.000) 


1.000(0.166) 


40 


1.000(0.000) 


1.650 


1.283(0.017) 


1.364(0.201) 


40 


1.200(0.043) 


2.200 


1.093(0.026) 


1.223(0.164) 


40 


0.940(0.038) 


2.469 


0.707(0.030) 


0.794(0.111) 


40 


0.551(0.076) 


2.510 


0.683(0.037) 


0.796(0.107) 


40 


0.493(0.050) 


2.550 


0.758(0.037) 


0.843(0.123) 


40 


0.511(0.040) 


2.591 


0.773(0.033) 


0.903(0.115) 


40 


0.547(0.002) 


2.631 


0.811(0.029) 


0.889(0.128) 


40 


0.518(0.041) 


2.671 


0.823(0.037) 


0.914(0.132) 


40 


0.532(0.032) 


2.712 


0.776(0.038) 


0.926(0.114) 


40 


0.492(0.035) 



a Thc complete vesrion of this table is in the electronic edition of 
the Journal. The printed edition contains only a sample. The average 
spectra are shown in /„, all normalized to have 1 Jy at J. The standard 
deviation of the mean is given in the parentheses following each flux. 

b Individual spectra normalized at J. 

individual spectra normalized by AFE(7.7). 

d Numbcr of galaxies used in the averaging. 

°The average spectrum of the 2 ellipticals using the J-band normal- 
ization scheme. Note that the AFE(7.7) normalization scheme would 
lead to essentially the same average spectrum for the ellipticals. 
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Table 4. Infrared Lines and Emission Features 



Wavelength 


Identifications 


EW 


(fims) 




H 


3.30 


AFE(3.3) 


0.02 


4.03 


Bra? 


0.06 


6.2 


AFE(6.2) 




7.0 


[Aril] 6.99/im + AFE(6.9) + H 2 0-0 S(5) 


0.003 


7.7 


AFE(7.6/7.7) 




8.6 


AFE(8.6) 




10.6 


[SIV] 10.51/im 


0.002 


11.3 


AFE(11.3) 





Table 5. Median NIR-to-FIR Flux Ratios 



Subsample 


No." 


R^O/lOO) 6 


log L Fm /L B b 


NIR total /FIR ft 


NIRjSM/ FIR 6 


NIR^M/ FIR 6 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


FIR-quiescent 


12 


0.39 (0.06) 


-0.18(0.24) 


0.174 (0.10) 


0.028 (0.03) 


0.041 (0.03) 


Intermediate 


11 


0.52 (0.08) 


0.37(0.35) 


0.083 (0.06) 


0.019(0.02) 


0.035 (0.02) 


FIR-active 


13 


0.73 (0.09) 


0.54(0.39) 


0.043 (0.04) 


0.013(0.03) 


0.020 (0.03) 



a Actual number of galaxies in each subsample. 

b Mcdian ratios, each followed in the parentheses by the r.m.s. dispersion with respect to the median. Columns (5) 
to (7) are the ratios of the integrated PHT-S flux over 2.5 to 4.7 (im to the FIR flux by treating the stellar continuum 
differently: no stellar continuum is removed in Column (5), and it is removed in both Columns (6) and (7) using 
the reddening corrected and zero-reddening methods, respectively. 



Table 6. Median Flux-Density Ratios of the AFEs 



No. a log L FIR /L b B (6.2)/(7.7) 6 (8.6)/(7.7) 6 (11.3)/(7.7) b 



(1) 


(2) 


(3) 


(4) 


(5) 


13 


-0.21(0.16) 


0.69(0.12) 


0.48(0.07) 


0.46(0.13) 


13 


0.20(0.15) 


0.69(0.05) 


0.48(0.08) 


0.47(0.11) 


13 


0.82(0.27) 


0.65(0.04) 


0.46(0.04) 


0.40(0.09) 


All 


0.20(0.46) 


0.66(0.09) 


0.47(0.06) 


0.45(0.12) 



a Number of galaxies in each subsample. The last row is for the whole 
sample. 

b Median ratios, each followed in parentheses by the r.m.s. dispersion 
with respect to that median. 
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Table 7. Median MIR-to-FIR Flux Ratios 



Subsamplc" 


(5.8-11.3)/FIR b 


AFEs(6-9)/FIR 6 


AFEs/FIR 6 


(1) 


(2) 


(3) 


(4) 


FIR-quiescent 


0.18(0.06) 


0.13(0.05) 


0.15(0.05) 


Intermediate 


0.17(0.04) 


0.13(0.04) 


0.15(0.04) 


FIR-activc 


0.11 (0.05) 


0.09 (0.04) 


0.09 (0.04) 



a Same subsamplcs as in Table 5. 

b Column (2) is the median ratio of the integrated flux between 5.8 and 
11.3 (im to the FIR flux. Column (3) is the median ratio of the summed 
flux of AFE(6.2), AFE(7.7) and AFE(8.6) to the FIR flux. Column (4) 
is the median ratio of the summed flux of all the 4 MIR features to the 
FIR flux. Each median flux ratio is followed in the parentheses by the 
r.m.s. dispersion with respect to that median. PHT-S aperture correction 
has been applied to all the MIR fluxes. 



Table CI. Feature Profile Slopes 



Profile Slope" 


Median 


Dispersion 6 




(1) 


(2) 


(3) 


(4) 


S- (6.2) 


0.69 


0.14 


0.05 


5+ (6.2) 


0.42 


0.14 


0.06 


S- (7.7) 


0.49 


0.08 


0.03 


5+ (7.7) 


0.33 


0.08 


0.03 


S- (8.6) 


0.09 


0.12 


0.04 


5+ (8.6) 


0.55 


0.15 


0.10 


S- (11.3) 


0.58 


0.24 


0.11 



a S'_ and S+ refer to the profile slopes at the short- 
arid long-wavelength sides of an emission feature, re- 
spectively. 

b The r.m.s dispersion with respect to the median. 

c Median of the statistical measurement errors in 
the profile slope. 
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Fig. I. — Distribution of the galaxies in a plot of logL F i R /L B vs. 72(60/100). The sample 
is divided into an "FIR-quiescent" subsample (open squares), an "FIR- active" subsample 
(solid squares), and an intermediate subsample (crosses). NGC 3620, for which Lb is not 
available, was assigned to the intermediate subsample, and is omitted from Figs. 1, 7, and 
8. 
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Fig. 2. — PHT-S spectra in units of Jy as a function of the rest-frame wavelength in microns. 
The error bars show the statistical errors for each spectral point. 
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Fig. 2. — Continuation of Fig. 2a. 
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Continuation of Fig. 2a. 
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Fig. 3. — Plots of the average rest-frame spectra derived from weighted averages of the 
spectra of 40 galaxies (see §3.1). The squares represent the average spectrum obtained 
by normalizing to the J fluxes, while the thick solid curve results from normalizing by the 
integrated flux of AFE(7.7 fim). Representative error bars are shown only for the former; note 
that the error bars should smallest near the fiducial wavelength used for the normalization. 
The spectra of two elliptical galaxies are also shown: NGC 3379 (the thin solid curve) and 
NGC4374 (the dotted line). 
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Fig. 4. — Comparison of the average spectra of the FIR-quiescent and FIR-active subsamples 
(see §2.1; also Table 5). The squares represent the FIR-active subsample and the solid curve 
is the average spectrum of the FIR-quiescent galaxies. The mean of the elliptical galaxy 
spectra is shown as a dotted curve. (Numerical values for all three curves can be found 
in Table 3.) As explained in the text, different normalizations are used for generating the 
spectra in (a) and (b). 
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Fig. 5.— Near-IR "color-color" plots: (a) log F X (K)/F X (J) vs. log F X (H)/F X (J), and (b) 
log F x (4/j,m)/F x (J) vs. log F X (H)/F X (J), where F A (4/im) is the mean flux density at 4/xm 
as defined in Table 2. The Key Project galaxies are shown as open squares, solid squares, 
and crosses (as in Fig. 1), and the reference elliptical galaxies as circles. In each panel, the 
solid line is the inferred reddening line for the case of a foreground dust screen the dotted 
curve is for the case where dust and stars are uniformly mixed. The tick marks along the 
solid line indicate respectively Ay = 1,2, and 3 m for the dust screen case. 
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Fig. 6. — Averaged spectra obtaining by subtracting from each disk galaxy spectrum (a) 
a stellar continuum reddened to agree with the H/J color, or (b) an unreddened stellar 
continuum (derived from the elliptical galaxies.) The solid curves are modified black-bodies 
with an emissivity that scales as A~ 2 and T = 750 K in (a) and 10 3 K in (6). The error bars 
represent the standard deviation of the mean. The individual, stellar continuum-subtracted 
spectra were normalized by the integrated flux of AFE(7.7) prior to averaging, and the results 
are renormalized to the same peak flux of the 7.7 /im feature. 
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Fig. 7. — Plot of the ratio of the flux density at 4/xm to the mean flux density AFE(7.7), as 
a function of Lfir/^b- Both F\(4/im) and (7.7) were computed after removal of the stellar 
contributions. The dotted line indicates the median sample flux ratio of 0.11. 
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Fig. 8. — Plots of the relative strengths of the AFEs as a function of Lfir/Lb, where 
(6.2)/(7.7), (8.6)/(7.7), (11.3) / (7.7) represent the mean Fx for these features after removal 
of the stellar continuum (using the zero- reddening method). We have divided the data set 
into three equal-size subsamples in terms of Lfir/Lb, and plotted their median values as 
large crosses in each plot. The extent of each cross indicate the r.m.s. dispersion with respect 
to the median value for each subsample. A few outliers are labelled. 
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Fig. 9. — Plot of the ratio of the combined luminosity of the AFEs to the FIR dust emission, 
as a function of i?(60/100). Galaxies with (without) a measured interstellar far-UV radiation 
field G m Malhotra et al. (2001) are represented by squares (circles). The sizes of the squares 
represent the approximate value of Go, and the filled symbols represent galaxies with the 
largest aperture corrections, as described in the text. The statistical errors in -^afes/^fir 
are small and thus no error bars are plotted. Note that nearly all of the galaxies fall within 
a wedge defined by the solid and dotted lines, derived from IRAS observations of reflection 
nebulosity and H II regions respectively. 



-44- 



1.5 



(X? 



CV! 



0.5 



□ □ 

□ n 



□ □ 



□ □ □ 



■■n% D 



0.2 0.4 0.6 0.8 1 

Fractional CAM LW2 flux within PHT-S aperture 



Fig. 10. — Comparison of the PHT-S flux scale with the flux scale of CAM in the LW2 filter 
(Ao = 6.7/im). The abscissa is p, the fraction of the galaxy flux at 6.7//m that falls within 
the PHT-S 24" x 24" aperture. The ordinate is the ratio of the CAM LW2 flux within the 
PHT-S aperture to a CAM LW2-band-equivalent flux derived from the PHT-S spectrum. 
The dotted line indicates a ratio of 0.76, the mean value for the 20 sources with p > 0.6. 
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Fig. 11. — Similar to Fig. 10, but comparing with CAM LW1 (4 to 5 /xm) data. The dotted 
line indicates a flux ratio of 0.7, the median value of the data shown in the plot. 
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Fig. 12. — Plots of (a) the flux density at 4/im and (b) the FIR-to-blue luminosity ratio as 
a function of K-band flux density for 40 galaxies. The galaxies with logLpm/^B < 0.1 are 
represented by solid squares; notice that these galaxies have a flatter distribution in both 
plots, especially in (b), where they show no trend with increasing F^. The dotted line in (a) 
is a least-squares fit to all 40 galaxies, by minimizing their distances perpendicular to the 
fit. The solid line is a similar fit to the galaxies indicated by solid squares. 



